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Equations (4) and (ll> are next substituted into 
equation (b). Further the convective tens of the 
substantial derivative in equation (b), ufos/cM 3 ), is 
neglected so that the final equations will have no 
nodes propagating at the flow velocity. Thus 
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where D/lk 3 is the substantial derivative D/D^ 

- + w 0 '/'X. 

Equations (l) to (?) are linearized by consid- 
ering small perturbations of the variables irons their 
equilibrium values. Substitution ot the s^all per- 
turbation form of the variables c, u, p, t, and s 
and e l iminati or. os squares and cross -products yields 
the following linearized equations written in non- 
dimensional form. 
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Consequently, the bulk gas temperature is related to 
the pressure by 



(13) 


Taking the time Fourier transform of equations (12) 
and (13) yields the following transfer functions re- 
lating a perturbation in air temperature and entropy 
to a perturbation in pressure: 
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Substituting equations (19) and (20) into equation (IS) 
yields the following Fourier transformed velocity po- 
tential wave equation: 



k s - u>/e s (23) 


and 

M * u /c (24) 

SOS 

The velocity potential solution is assumed to be pro- 
portional to exp f"(ik ii) - ia*M. Substituting this 
solution into equation (21) yields a wave number equa- 
tion which has as a solution: 
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APPLICATION TO BOCTED COWS TIC* STSTOf 

The analysis developed in the last section is 
now used to study the spectral structure of pressure 
measurements made ia the ducted combustion system 
shown in Fig. i. The ducted combustion system shown 
in Fig. I consists of: (l) a source region inside 

the combustor can, ( 2 ) a non-source region inside a 
spool piece and a long duct, (3) an area expansion 
and contraction on either side of the spool piece, 
and (4) an upstream boundary at the exit of the long 
duct and a downstream boundary at the combustor inlet. 

The solution for the velocity potential in a non- 
source region is given by equation (21). Th Fourier 
transformed velocity potential wave equation (eq. (21)) 
is assumed to apply in the source region with the 
addition of a source term C(a,x) on the right-hand 
side. I'sing the coobustor inlet and exit impedance 
as a boundary condition, a unique Green’s function 
solution for the Fourier transformed velocity poten- 
tial is found. The Fourier transformed acoustic 
pressure and particle velocity can be found from a 
velocity potential using equations (19) and (20). 
Consequently, the acoustic pressure and particle vel- 
ocity at the combustor exit can be determined from 
the velocity potential solution in the combustor. The 
acoustic pressure 3nd particle velocity at the com- 
bustor exit can then be used to find the acoustic 
pressure and particle velocity at any other point in 
the ducted combustor system with four -pole equations 
in transfer matrix form derived using equations (19), 
(20), and (27). The combustor exit impedance can be 
calculated from the duct exit impedance using the same 
four-pole transfer matrix equations. The source re- 
gion solutior, duct transfer matrix, area discontin- 
uity matrix, and boundary conditions are discussed in 
more detail in reference (1). 

The following procedure is used to calculate the 
pressure spectrum. It is applied at each frequency 
as necessary to obtain the desired spectrum. First, 
the exit pressure perturbation is arbitrarily assumed 
to be one Pascal. Then using the duct exit acoustic 
impedance the particle velocity is calculated. Next, 
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and the cross -spectrum between tun points is calcu- 
lated trom 
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ItK eiperUMtiul «pp«ratus Is sboiasdMHCicallj 
ia Fig. 1. combustor section consist s of a J-47 

burner cm placed concentric* I ly in a 0.30m diameter 
by 0o?7n long flow duct. The combustor section is 
follwed by a 0.38* diameter *>y 0.76a long spool 
piece. This section is followed by a 0.30m diao- 
eter by 0.1m long flow duct. 

The measurements discussed herein were made at 
an exit temperature of 920 K and at air mass flow 
rates of 0.5, 1.13, and 1.6$ kg/sec. The correspond- 
ing velocities at the exit of the long duct were 
18.5, 41.6, and 61.3 m/sec and the corresponding fuel 
flow rates were 0.009, 0.01$, and 0.027 kg/sec. The 
fuel used was Jet A. 

Simultaneous internal fluctuating pressure mea- 
surements were made at the three locations shown in 
Fig. I. The transducers used ere conventional 5/8 
cm diameter (nominal) pressure response condenser mi- 
crophones. To avoid direct exposure to the severe 
environment within the flow duct, the microphones 
were mounted outside the duct and the fluctuating 
pressure in the duct was communicated to the trans- 
ducers by ' ’semi-infinite' * acoustic waveguides. 

More information on the probes is given in refer- 
ence (2), 

Constant-bandwidth pressure auto-spectra meas- 
ured at three test conditions are shown in Fig. 2. 

The spectra measured in the spool piece upstream of 
the area contraction are shown in Fig. 2(a). The 
spectra measured near the entrance to the long duct 
downstream of the area contraction are shown in Fig. 
2(b). The structure of the spectra shown in Figs. 
2(a) and 2(b) does not change greatly with operating 
condition over the frequency range from 2 to 200 Hz. 
Also the location of resonance peaks and dips in 
Figs. 2(a) and 2(b) is nearly the same for each oper- 
ating condition. Moreover, the peaks tend to be 
sharper at the low frequencies and more broad at the 
higher frequencies. However, the spectra shown in 
Fig. 2(a) have a different structure than those shown 
in Fig. 2(b). 


Cams t amt hamArtdtb pressure cross -spectra meas- 
ured across tbe area com tract tom at three test condi- 
tion* are th orn s, in Figs. 3(a), 3(b), and 3(c), re- 
spectively. Both nmgnltamde and phase are shorn. The 
location of the r esonance peaks and dips in the mag- 
nitude pint of the crons spectra is nearly ~he sane 
at each test condition. Also, the peaks teed to he 
sharper at the law frequencies and more broad at the 
higher frequencies. The phase plat c —at be char- 
acterized as either Just a negative sloped straight 
lime due to a tine delay nor just a curve chat varies 
abruptly between 0 and ISO degrees ia a regular man- 
ner that would be due to the presence of wares having 
appro ximately equal magnitude but opposite sign tra- 
velling upstream and downstream. Instead the phase 
plot shew; both tendencies with the tine delay ten- 
dency dominating at the higher frequencies. 

A large peak appears above 200 Hz in die spectra 
measured at 41.6 a /see and at 61.3 m/sec. These 
peats appear at such high frequencies that they can- 
not be related to specific duct resonance harmonics. 
The pedes nay be due to a feedback resonance between 
the combustor and the duct. They are not analyzed in 
this paper. 

COKPARISOK WITH EXPEKB&HTAL DA I A 

Results determined by employing the model using 

T - 1.0 and x /* ■ 0.001 are shown in Figs. 

s s s s 

4 to 7. Exhibited in Fig. 4 is the sound propagation 

velocity and the attenuation. For T /< =1-0 the 

s s 

attenuation is zero and the sound propagation veloc- 
ity is equal to the adiabatic speed of sound. How- 
ever, for t /*• - 0.001 the attenuation is greater 

s s 

than 0.5 dB from 50 to 320 Hz and the sound propaga- 
tion velocity varies from 50$ to 580 m/sec. 

The auto-spectra and cross -spectra calculated 

for t A «• 1.0 and x /< - 0.001 are compared 

s s s s 

with the data measured at v ■ 18,5 a/sec in Figs. 5 

to 7 . The model calculations made using x /* 

s s 

= 1.0 produce results which do not agree with the 
data. However, the calculated auto-spectra appear to 
be stretched out versions of the measured auto- 
spcctra. 

The model calculations of pressure auto-spectra 
and cross-spectra made using x^/* s * 0.001 show bet- 

ber agreement with the measured data. The peaks and 
dips in the calculated pressure auto-spectra and 
cross-spectra are closer to the locations of the peaks 
and dips in the measured data. In addition, the 
structure of the measured and calculated curves is 
similar. 

DISCUSSION 

The large discrepancy between the auto-spectra 

calculated using the model with x /< *1.0 (no at- 

s s 

tenuation and a sound propagation velocity of 610 
m/s) and the measured pressure auto- spectra compared 
in Figs. 5(a) and 6(a) might have two explanations . 
First, the approach used to model the system might 
not apply. For example, the dips and peaks in the 
measured spectra could be due to some complicated 
feedback interaction process between the combustion 
noise source and the duct which would not be related 
to the manner in which sound propagates in the duct. 
Second, the approach could be correct and the values 
of attenuation and sound propagation velocity mi ght be 
in error. The results obtained using the model with 
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the * 0-001 sknn is Figs. S(k) and b(b) wj- 

g«t chs scssscic wdcl is correct sines Che charac- 
teristic structure of the measured spectra between 2 
and 200 Hz is similar to the calculated spectra. 
However, the best argument tV; the correctness of the 
acoustic node l is the agreement between the measured 
and calculated yfcase angle curves compared in Fig. 
7(h) since the cross -spec true phase angie is indepen- 
dent of the source region. While these results indi- 
cate that the acoustic node l is correct, they also 
suggest that a more coup lex model for the soot parti- 
cle oxidation process nay be necessary, the use of a 
constant time- independent soot particle temperature 
greatly over simplifies a complex problem, in addi- 
tion, the combustor inlet and duct exit boundary 
conditions used herein may need to be improved. 

C0NCITDIHC REMARKS 

The m>del discussed provides a good method for 
predicting pressure auto-spectra and cross-spectra 
at an area contraction. The method depends only on 
a single parameter. 

The results obtained indt -ate that dispersion 
and attenuation may be a fairly cc ma i tm occurrence in 
liquid fuel combustor test facilities., m* i v suggest 
that it may not be correct to assume th** *ovmd pro- 
pagates adiabaticaily in combustor test tacilities. 

The propagation of sound in a combustion duct 
facility appears to be a complex process. The pre- 
sence of oxidizing soot particles may be affecting 
slcongly the propagation of sound at combustion noise 
frequencies. Theoretical analysis of the affect of 
oxidizing soot particles is in a primitive stage of 
development. Further theoretical analysis seems 
warranted. Also, experiments which are designed to 
ascertain functional dependences that can be checked 
against theory and which include both acoustic and 
coobustion measurements would be of interest . 
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Figure !. Ducted combustion system. 
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(a) MEASURED IN SPOOL PIECE UPSTREAM OF AREA CONTRACTION. 



(b) MEASURED NEAR ENTRANCE TO LONG DUCT DOWNSTREAM 
OF AREA CONTRACTION. 


Figure 2. - Pressure auto-spectra. 
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(c) V *61.3 m/sec. 


Figure 3. - Measured pressure cross-spectra across area 
contraction. 
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(b) Tj/k j ■ 0.001. 


Figure 5. - Comparison of measured and calculated auto-spectra 
in spool piece downstream of area contraction (V • 18.5 m/sec). 
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Figure 6. - Comparison of measured and calculated auto-spedra 
near long duct entrance upstream of area contraction (V * 18. 5 
m/sec). 
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